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Abstract 



In a recent paper (Phys. Rev. B 59, 9699 (1999)), Chandross and Hicks 
claim to present a new density matrix renormalisation group (DMRG) method 
for dealing with excited states of quantum lattice models. The proposed im- 
provement to the DMRG — the inclusion of excited state wave functions in ad- 
dition to the ground state in the density matrix when calculating excitations — 
is in fact standard pratice, is clearly stated in White's original papers, and has 
been used repeatedly by many groups to study excited states. The authors 
apply the method to the extended, dimerised Hubbard model for conjugated 
polymers. The criteria for determining whether states are bound or not are 
assessed. The authors claim that their results show that the optically im- 
portant "l-B^i" state is bound (excitonic), in contrast to a previous study. 
However, the discussion is qualitative, and the authors arrive at conclusions 
on the basis of results for one lattice size only. We show that when Chandross 
and Hicks' criterion is developed into a quantitative definition of particle-hole 
separation, with the finite-size dependence analysed, the implication is that 

the \Bu state is unbound, in keeping with the conclusions of a previous study. 
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In a recent paper [|l| , Chandross and Hicks claim to present a new density matrix renor- 
malisation group (DMRG) method [0,^ for dealing with excited states of quantum lattice 
models. They apply the method to the dimerised, extended Hubbard model for conjugated 
polymers. They claim that a previous study P of this model is flawed because it uses a 
"conventional" DMRG method which does not handle excitations correctly. The improve- 
ment that they suggest is to form a density matrix not only from the ground state, but 
from all the states being targeted in the calculation. This is in fact standard practice in 
DMRG calculations of excited states and the structure of the density matrix required to 
target excited states is given in White's original papers on the method It has been 

used by many authors to target excitations in a variety of quantum lattice models (see, e.g., 
0) and was certainly used in when various excitation energies and correlation functions 
were calculated for the extended Hubbard model. The comparisons presented in Fig. 1 and 
Fig. 2 of between the "conventional" DMRG and Chandross and Hick's "improvement" 
are therefore of limited value, as, to the best of our knowledge, all DMRG studies of excited 
states to date have incorporated the targeted excitations into the density matrix [|^ . Unfor- 
tunately, a slightly different value for the Coulomb V is used in so a direct comparison 
with the results (e.g., for energies) tabulated in is not possible. We have run a DMRG 
program which uses the algorithm used in |^ for targeting excited states with the param- 
eters U = 3t, V = t, 6 = 0.1, used in [|l|, and found good agreement for the energies and 
correlation functions with the results plotted in Fig. 1(a) and Fig. 2(a) of [|l|. For instance, 
we plot the IB^ and niAg energies as functions of the lattice size in Fig. |l|. The results 
compare very well with Fig. 1(a) of |jl[. 
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FIGURES 




N 

FIG. 1. The energies of the IB^ (diamonds) and mAg (triangles) states of the dimerised, 
extended Hubbard model as a function of the lattice size N for the parameter set used in Q . The 
number of states retained per block [0,0] is m = 270. 



In []T[] Chandross and Hicks also examine criteria for deciding whether a particular ex- 
citation is bound (excitonic) or not. They claim that the average particle-hole separation, 
defined in in terms of the density-density correlation function, is "too approximate" a 
quantity to determine whether a state is bound or not. They argue that by inspecting the 
centered correlation function as a function of distance (together with the profile of doubly 
occupied sites along the chain), for one particular lattice size (A^ = 36 sites), one can see that 
the IBu and mAg states are "different" in that the IBu {mAg) has its strongest particle-hole 
correlations at short (long) distances. However, Chandross and Hicks do not present an 
alternative quantitative definition of the particle-hole separation, based on this observation. 
In , on the other hand, it is argued that a definition of particle- hole binding must take into 
account the way in which correlations scale with lattice size A^. In it is argued that this 
scaling is different for bound and unbound excitations, and that the scaling of the average 
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particle-hole separation with is but one manifestation of this. 
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FIG. 2. The average number of doubly occupied sites of the IB^ state relative to the ground 
state at distance i from the center of the chain for various lattice sizes N. 

Suppose we wish to take the average double occupancy of the IB^ state (relative to the 
ground state) along the chain {{hi — 1)^)ib^ — {{hi — l)^)iyig as an example (Fig. 2(a) in 
||l|]). In Fig. H we plot this quantity for various lattice sizes A^. We see that, although the 
concentration of doubly occupied sites is greatest in the middle of the chain, the distribution 
spreads out as is increased. The area under these curves rapidly converges to a non-zero 
value (~ 0.538) as ^ oo. This shows that the number of pairs of particles and holes in 
the IBu, relative to the number in the ground state, approaches a constant. Our results 
could indicate that particle-hole pairs separate as A^ is increased and are hence unbound, or 
they may simply indicate dispersion of a bound exciton in the IB^- 
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FIG. 3. The averaged, centered, odd-site correlation function (relative to the ground state 



(solid diamonds). 

To address this, we again consider the averaged, centered, odd-site correlation function 
C^^(z), (again relative to the ground state value), defined in |^ and plotted for = 36 in 
the inset to Fig. 2(a) in In Fig. ^we plot this quantity for a number of values of A^. We 
see that, although the correlations are generally strongest at short distances, they become 
increasingly spread out, and hence the particle-hole pair becomes increasingly separated, as 
is increased. Indeed, if one utilises |C'j^^(j)| to define a probability distribution for the 
particle- hole separation, as in 0, then one finds that the average particle- hole separation 
grows linearly with A^, as shown in Fig. ^. We note that any use of the density- density 
correlation function to describe particle-hole separation and the nature of exciton binding 
of excited states in the extended Hubbard model is merely plausible rather than rigorous 
0], but Chandross and Hicks [|T| do not offer an alternative quantitative definition of the 
particle- hole separation to the ones provided in 



value, as defined in ||5|) for the IB^ state for N = A2 (diamonds), 62 (triangles), 82 (stars) and 102 
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FIG. 4. The (reduced) average particle-hole separation, as defined in |5[ by using |C^^(i)| as 
a probability distribution, for the IBu state. Note the linear increase with A'^. 

To summarise, Chandross and Hicks claim that, because the IB^ and mAg have their 
greatest particle-hole correlations at short and long distances respectively (on the = 36 
lattice), the IB^ is bound and the mAg is unbound. We would argue that it indicates that 
the particle-hole potential is more strongly attractive for the state than for the mAg. 
However, from the plausible, quantitative definition of the particle-hole separation given 
above, it would appear that the attraction between the particle and hole in the state 
is not strong enough to bind them, and their separation increases throughout the range of 
lattice sizes studied. 

Finally, we consider the structure of the density matrix when targeting excitations such as 
the mAg and nBu- Chandross and Hicks argue that only four states — the lAg (ground state), 
the IBu, the mAg and the nB^ — need be included in the density matrix. Our examinations 
of the dipole moments between the Ag states and the IB^ indicate that this approach is 
probably reasonable for the mAg which is well defined. That is, there is a reasonably abrupt 
jump in the magnitude of the dipole moment {lBu\fi\jAg) at j = m. As shown in [||, this 
coincides with jump in the ionicity (the average number of doubly occupied sites) and in 



the particle-hole separation. However, the nBu state is less well defined in that there can be 
a number of excitations that have a strong dipole moment with the mAg. This can be 
seen in Table | where we list the dipole moments {jBu\fl\mAg) for = 6, 10, 14 and 18, for 
the first five Bu states. Note that in no case is the nB^ state clearly defined, though there 
is a general trend whereby the 25„ increases its relative dipole strength with the mAg at 
the expense of the 4S„. Our contention here, as proposed in is that, at least in terms of 
dipole moments or the density-density correlation function, the IBu state is the threshold 
of unbound states in the Bu sector and the ^^nBu" is not well defined for this model. 

Calculations were performed at the New South Wales Center for Parallel Computing. 
This work was supported by the Australian Research Council. 



7 



REFERENCES 



* Email address: phlrb@newt.phys.unsw.edu.au 

[1] M. Chandross and J. C. Hicks, Phys. Rev. B 59, 9699 (1999). 

[2] S. R. White, Phys. Rev. Lett. 69, 2863 (1992). 

[3] S. R. White, Phys. Rev. B 48, 10345 (1993). 

[4] Note also that the dipole operator method for calculating dip ole- allowed singlet states 
within the DMRG proposed in |l[|, was actually first described in G. P. Zhang, T. F. 
George and L. N. Pandey, J. Chem. Phys. 109, 2562 (1998). 

[5] M. Boman and R. J. Bursill, Phys. Rev. B 57, 15167 (1998). 

[6] S. R. White and D. A. Huse, Phys. Rev. B 48, 3844 (1993). 

[7] We note further that Chandross and Hicks |T]] claim that including only the ground 
state in the density matrix leads to "grossly incorrect quantitative" results for excited 
states. Although it is illogical to include only the ground state in the density matrix when 
calculating excited states, for a given lattice size (e.g. = 36) it is in fact a mathematical 
certainty that the DMRG results for the excitations with a simple density matrix formed 
solely from the ground state projection operator will converge to the exact results as the 
number of states retained per block, m, is increased. The phrase "less rapidly converging" 
is therefore more appropriate, as DMRG results should not be accepted if they have not 
converged with m. 

[8] The mAg state is defined to be the excited state in the ground state symmetry sector 
that has the largest dipole moment with the IBu state, which is the lowest energy state 
in the dipole allowed symmetry sector. 

[9] M. Yu. Lavrentiev and W. Barford, Phys. Rev. B 59, 15048 (1999). In this paper an alter- 
native quantitative definition of the particle-hole separation is proposed using molecular 



8 



orbital operators. In terms of this definition there are indications that the 1S„ particle- 
hole separation grows less rapidly than linearly with N. 



9 



TABLES 







i = l 


i = 2 


i = 3 


j = 4 


j = 5 


6 




2.32 


0.76 


0.30 


1.87 


1.39 


10 




3.48 


1.77 


0.38 


3.30 


0.06 


14 




4.45 


3.15 


0.04 


3.98 


0.10 


18 




5.33 


4.73 


0.67 


4.24 


1.73 




TABLE 1. 


Transition 


moments with the 


mAg states 


for the first five 


states (i.e. 



{j By,\jl\mAg) for j = 1,. . . ,5) for AT = 6, 10, 14 and 18. Note that there is no clearly defined 
"nB" state. 
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